Specifics of impurity effects in ferropnictide superconductors 
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Effects of impurities and disorder on quasiparticle spectrum in superconducting iron pnictides 
are considered. Possibility for occurrence of localized energy levels due to impurities within the 
superconducting gap and the related modification of band structure and of superconducting order 
parameter are discussed. The evolution of superconducting state with impurity doping is traced. 
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I. INTRODUCTION 

The recent discovery of superconductivity (SC) with 
rather high critical temperature in the family of doped 
ferropnictide compounds [1, 2], has motivated a great in- 
terest to these materials. Unlike the extensively studied 
cuprate family [3], that present insulating properties in 
their initial undoped state, the undoped LaOFeAs com- 
pound is a semimetal. The scanning tunnel microscopy 
(STM) study [4] established that this material has a lay- 
ered structure, where the relevant for SC layer is FeAs 
with a 2D square lattice by Fe atoms and with As atoms 
located out of plane, above or below the centers of square 
cells (Fig. 1). Its electronic structure, relevant for con- 
structing microscopic SC models, have been explored 
with high-resolution angle-resolved photoemission spec- 
troscopy (ARPES) techniques [5, 6]. Their results indi- 
cate the multiple connected structure of Fermi surface, 
consisting of electron and hole pockets and absence of 
nodes in both electron and hole gaps [5] , suggesting these 
systems to display the so-called extended s-wave SC or- 
der, changing its sign between electron and hole segments 
[10]. 

To study the band structure, the first principles nu- 
meric calculations are commonly used, outlining the im- 
portance of Fe atomic c?-orbitals. The calculations show 
that SC in these materials is associated with Fe atoms 
in the layer plane, represented in Fig. 1 by their orbitals 
and the related hopping amplitudes. The dominance of 
Fe atomic 3d orbitals in the density of states of LaOFeAs 
compound near its Fermi surface was demonstrated by 
the local density approximation (LDA) calculations [7- 
12]. It was then concluded that the multi-orbital effects 
are important for electronic excitation spectrum in the 
SC state, causing formation of two gaps: by electron 
and hole pockets at the Fermi surface. To explain the 
observed SC properties, it is suggested that these mate- 
rials may reveal an unconventional pairing mechanism, 
beyond the common electron-phonon scheme ]13, 14]. In 
general, the total of 5 atomic orbitals for each iron in the 
LaOFeAs compound can be involved, however the ways 




FIG. 1: Schematics of a FeAs layer in the LaoFeAs compound 
with dxz (white) and dyz (dark) Fe orbitals and the Fe-Fe 
hopping parameters in the minimal coupling model. Note that 
the hoppings between next near neighbors (is, 4) are mediated 
by the As orbitals (out of Fe plane). 



to reduce this basis are sought, in order to simplify an- 
alytical and computational work. Some authors ]15, 16] 
have suggested that it is sufficient to consider only the 
dxz and dyz orbitals. Thus, building such minimal cou- 
pling model based on two orbitals, one is able to adjust 
the model parameters (energy hopping and chemical po- 
tential) to obtain the Fermi surface with the same topol- 
ogy that found in the first principles calculations of band 
structure. 

Having established the SC state parameters, an im- 
portant class of problems can be considered about the 
effects of disorder, in particular by impurities, on the 
system electronic properties, and this issue has been also 
studied for doped ferropnictides. Alike the situation in 
doped perovskite cuprates, here impurity centers can ei- 
ther result from the dopants, necessary to form the very 
superconducting state, or from foreign atoms and other 
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local defects in the crystalline structure. Within the min- 
imal coupling model, an interesting possibility for local- 
ized impurity levels to appear within SC gaps in doped 
LaOFeAs was indicated, even for the simplest, so-called 
isotopic (or non-magnetic) type of impurity perturba- 
tion [17, 18]. This finding marks an essential difference 
from the traditional SC systems with s-wave gap on a 
single-connected Fermi surface, were such perturbations 
are known not to produce localized impurity states and 
thus to have no sizeable effect on SC transition tempera- 
ture, accordingly to the Anderson theorem [19]. In pres- 
ence of localized quasiparticle states by isolated impu- 
rity centers, the next important issue is the possibility 
for collective behavior of such states at high enough im- 
purity concentrations. This possibility was studied long 
ago for electronic quasiparticles in doped semiconducting 
systems ]20] and also for other types of quasiparticles in 
pnononic, magnonic, excitonic, etc. spectra under impu- 
rities ]21], establishing conditions for collective (including 
coherent) behavior of impurity excitations with striking 
effects in observable properties of such systems. As to 
the high-Tc doped cuprates, it is known that their rf-wave 
symmetry of SC order only permits existence of impurity 
resonances [22, 23], not the true localization, and hinders 
notable collective effects on their observable properties. 
As to our knowledgement, no consistent study on collec- 
tive impurity effects is know for the doped ferropnictide 
systems up to the moment, and this defines the main 
emphasis of the present work. Namely, we shall develop 
an analysis of these systems, using the Green function 
(GF) techniques, similar to those for doped cuprate SC 
systems ]24] , the minimal coupling model by two orbitals 
for ferropnictide electronic structure, and the simplest 
isotopic type for impurity perturbation. The structure 
of quasiparticle spectrum near in-gap impurity levels at 
finite impurity concentrations, conditions for emergence 
of specific branches of collective excitations in this region 
of the spectrum, and expected observable effects of such 
spectrum restructuring will be discussed. 



II. MODEL HAMILTONIAN AND GREEN 
FUNCTIONS 
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FIG. 2: Electron (— ) and hole (-)-) segments of the Fermi sur- 
face in the normal state of model system with electronic spec- 
trum by Eq. 5. The center of first Brillouin zone is displaced 
by (7r/2a, 7r/2a) to fully include all the segments around four 
characteristic points F, X, M, and Y in this zone. 



vectors 5x.y point to its nearest neighbors in the square 
lattice. Passing to the operators of orbital plane waves 
Xk,o- = Af~^^^ X^n^*"^ "^n-cr (with the number TV of lattice 
cells) and analogous ?/k,cr, and defining an "orbital" 2- 
spinor ?A^(k, cr) — (a;k,(T, yk.cr), one can expand the spinor 
Hamiltonian in quasimomentum: 



Ht 



k,(T 



(2) 



Here the 2x2 matrix 

/lt(k) = £ + ,kO-0 + £-M^i + SxyM^l (3) 

includes the Pauli matrices cTi and the energy functions 



with 



For the minimal coupling model of Fig. 1 , the hopping 
Hamiltonian Ht is written in the local orbital basis as: 



Ht 



E 



+ ^2 {x1^aXn+5y,a + yLayn+S^,a + h.C.) 

+ ts (s^n.o-S^n+^jc+^y.cr + xl^.cr'^n+S^ -Sy ,a 

+ Vn,aVn+S^+Sy,a + vLaVn+S^ -Sy .a + h.C.) 

+ ti (xl^aVn+S^+Sy.a + vLa^n+S^+Sy ,a 

- xiayn+S^-Sy,a - yi,aXn+S^ -Sy ,a + h.C.)] . (1) 

where Xn,^ and ?/„ are the Fermi operators for dxz and 
dyz Fe orbitals with spin cr on n lattice site and the 



£x,k = — 2^1 cos /cj: — 2^2 cos /cj^ — 4i3 cos /ca; cos /Cj^, 
Ej^^k — ~2ti cos ky — 2t2 cos — 4i3 cos kx cos ky, 
£xyM = ~ At i sin kx Sin ky. (4) 



An optimum fit for the calculated band structure 
within the minimum coupling model is attained with 
the following set of hopping parameters (in |ii| units): 
ti = —1.0, t2 = 1.3, ^3 = ^4 = —0.85, and with the choice 
of the Fermi energy (chemical potential at zero temper- 
ature) ef — 1-45 ]12]. The ht matrix is diagonalized by 
the standard unitary transformation: 



;7(k) 



_ ( cos6lk/2 -sin6lk/2 
sin6'k/2 cos6'k/2 



3 



with 0]^ = arctan (£a;j, k/£-,k)j transforming from the or- 
bital to subband basis: 



The energy eigenvalues in Eq. 4: 



£e,k 
Shk 



efc,e(k)=£+,k±-^£^^,k + £? 



>k' 



(5) 



(6) 



correspond to the two subbands in the normal state spec- 
trum that respectively define electron and hole pockets 
of the Fermi surface. There are two segments of each 
type, defined by the equations £e,/i(k) = /x, as shown in 
Fig. 2. We note that both functions cos6'k and sin6'k 
change their sign around these segments, corresponding 
to their "azimuthal dependencies" around characteristic 
points of the Brillouin zone (Fig. 2), so that integrals 
of these functions with some azimuthal-independent fac- 
tors over the relevant vicinity of Fermi surface practically 
vanish and are neglected beside such integrals of fully 
azimuthal-independent functions in the analysis below. 

The adequate basis for constructing the SC state is 
generated by the operators of electron and hole subbands: 

Q!k,<7 = Xk,^ cos0k/2 - yk,crsin6'k/2, 

/5k,cr = ?/k,(TCos6'k/2 + a:;k,<Tsin6'k/2, (7) 

giving rise to the "multiband-Nambu" 4-spinors "pj^ = 

^a][^ .|,,a_k,4,,/3k,t)/3-k,4,^ and to a 4x4 extension of the 
Hamiltonian Eq. 2 in the form: 



k,CT 



(8) 



where the 4x4 matrix 



/is(k) = /i(,(k) ig) fs + AkO-Q ig) fi, 

includes the Pauli matrices fi acting on the Nambu 
(particle-antiparticle) indices in \E'-spinors and hb{k) is 
defined by Eq. 5. The simplified form for the extended 

s-wave SC order is realized with the definition of the gap 
function by constant values, Ak = A on the electron 
segments and Ak = —A on the hole segments. 

The electronic dynamics of this system is determined 
by the (Fourier transformed) GF 4x4 matrices [21, 24, 
25]: 



the explicit GF for the unperturbed SC system with the 
Hamiltonian H^, Eq. 7, is diagonal in quasimomentum, 
Gk,k' = Su,k'Gl_ and 



Gl = 



+ 



efp +£e(k)f3 + Afi 

£To + £/t(k)f3 - Afi 



)cr+ 



(11) 



where a± = (ctq ±(73) /2 and the secular denominators 
-Di,k = e^— (k) — A^ for i = e,h. In what follows, we use 
the energy reference to the Fermi level £f and approxi- 
mate the segments of Fermi surface by some circles of 
radius ki around the characteristic points in the Bril- 
louin zone, so that the dispersion laws Sj (k) = ep + ^j,k 
permit to linearize the quasiparticle dispersion close to 
the Fermi level as ^j.k ~ f^'Vj (|k — Kj| — ki). Generally, 
the Fermi wavenumbers kj and related Fermi velocities 
Vj for j = e and h can somewhat differ at a given choice 
of hopping parameters and chemical potential, but, for 
simplicity, we shall neglect this difference and consider 
their single values kj = kp and vj = vp. 



III. IMPURITY PERTURBATION AND 
SELF-ENERGY 

We pass to the impurity problem where the above 
Hamiltonian is added by the perturbation terms due to 
non-magnetic impurities [17] on random sites p in Fe 
square lattice with an on-site energy shift V: 



Without loss of generality, the parameter V can be taken 
positive, and for GF calculations, this perturbation is 
suitably expressed in the multiband-Nambu basis: 

Himp = ^ E e*^'''-'') P*|,14,k'*k'. (13) 

p,k,k' 

through the 4x4 scattering matrix Vk.k' ~ VUlU]^' ® T3. 
Within the approach of Refs. [21, 24[, the solution for 
Eq. 9 with the perturbed Hamiltonian + Hi can be 
obtained in different forms, suitable for different types of 
states, band-like (extended) or localized. All these forms 
result from the basic equation of motion: 



(9) 

whose energy argument e is understood as £ — iO and 
{{A{t)^B {{))}) is the quantum statistical average with 
Hamiltonian H of the anticommutator of Heisenberg op- 
erators. From the equation of motion: 



£Gk,k' = Mk.k'O-o ®To+ (([*k,-H"] |*k'))' 



(10) 



p,k" 

by specific routines of iterating this equation for the 
"scattered" GF's Gk",k'. 
Thus, the algorithm, where the next iteration step 

never applies to the scattered GF's already present af- 
ter previous steps, e.g. that with k" = k in Eq. 14, 
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leads to the so-called fully renormalized form, suitable 
for band-like states: 



(G^)"'-Sk 



(15) 



where the self-energy matrix Sk is expressed by the re- 
lated group expansion (GE): 



Sk = cfk(l + cBk + ...) . 



(16) 



Here c = Yl^N ^ is the impurity concentration (per 
Fc site) and the T-matrix results from all the multiple 
scatterings by a single impurity: 



Tic = Vk,k + ^ t4,k'Gk/Vk',k 



IV. T-MATRIX AND QUASIPARTICLE STATES 

The T-matrix, Eq. 16, is readily simplified taking into 
account that 14;. k = ^o"o ® T3 and introducing the inte- 
grated Green function matrix: 

k 

This diagonal form follows directly from the aforemen- 
tioned cancellation of the integrals with cos 6^ and sin ^k 
that appear in all the matrix elements of Ui^G^U^ except 
those proportional to S(t± (g) fo- Respectively, the func- 
tions gj{e) = N^^ J2k ^Jk i = e, /i are approximated 
near the Fermi level, |£ — ep] ^ as: 



k'^^k 



Tip, 



(20) 



-|- 14,k'Gk't4;',k''Gk''t4'',k + • (17) Here = TOja^/(27rfi,^) are the Fermi densities of states 



Ar2 

k'#k,k'Vk,k' 

The next term to the unity in the brackets in Eq. 14: 

Sk = X (i„e-* " + i„i-„) (1 - i„i-n)"' , (18) 
11 

describes the effects of indirect interactions in pairs of 
impurities, separated by vector n, in terms of interaction 
matrices A„ = ^k X^k'^^k ^''^ "Gk'- Besides this restric- 
tion on summation, multiple sums in the products like 
^n^-n never contain coincident quasimomenta. Eq. 18 
presents the first non-trivial GE term and the rest of its 
terms omitted in Eq. 14 correspond to the contributions 
from groups of three and more impurities [21]. 

An alternative iteration routine for Eq. 13 applies it 
to all the scattered GF's, this results in the so-called 
non-renormalized form, suitable for localized states: 



for respective subbands (in parabolic approximation for 
their dispersion laws), and by the assumed identity of all 
the segments of Fermi surface they can be also considered 
identical pj — pp. Omitted terms in Eq. 16 are of higher 
orders in the small parameter |£|/£f ^ 1- 

Then the momentum independent T-matrix is explic- 
itly written as 



7 



eoT3 



=■2 _ =-2 



(21) 



Gk 



AO fiO AO 
'^k^k'^k- 



(19) 



Here the non-renormalized self-energy GE: = 
cT ^1 + C-B^ -f . . . ^ , differs from the above renormalized 
one by absence of restrictions in quasimomentum sums 
for interaction matrices = Tk J2k' "G^' and their 
products. 

At the first step, we shall restrict GE to the common T- 
matrix level, providing the conditions for localized quasi- 
particle states with in-gap energy levels to appear at sin- 
gle impurities [16], and will study certain (narrow) energy 
bands of specific collective states that can be formed near 
these levels at finite impurity concentrations. At the next 
step, the criteria for such collective states to really exist 
in the disordered SC system will follow from the analy- 
sis of non-trivial GE terms. We notice that presence of 
renormalized GF's Gk' in the above interaction matrices 
is just necessary for adequate treatment of interaction 
effects over the in-gap bands. 



where £0 = A/vT+i? defines the in-gap impurity level 
[16] through the dimensionless impurity perturbation pa- 
rameter V = irppV, and 7^ = v'^Veq/A is the effective 
constant of coupling between localized and band quasi- 
particles. 

At finite c, using this T-matrix in Eq. 14, we ob- 
tain, from the condition detGj^"^ = [25], the formal 
dispersion equation expressed through dispersion of nor- 
mal quasiparticles ^k = £k — £f (but in neglect of the 
energy level width due to the effects of indirect interac- 
tion between impurities by higher GE terms): 



0. (22) 



Its solutions shown in Fig. 3 in function of the quasimo- 
mentum argument ^ = ^k display a peculiar multiband 
structure. First of all, it includes four modified bands 
±£f,(±^), slightly shifted with respect to the unperturbed 
SC quasiparticle bands i-^/A^ -|- accordingly to the 
basic function: 



£6(0 « V^Te 



■ C7 



V^+eie+eo)' 



(23) 



with = A^ — £q. It should be noted that these sub- 
bands for opposite signs of their argument ^ in fact refer 
to excitations around different segments (by electron and 
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FIG. 3: Dispersion laws for band-like quasiparticles in the 
T-matrix approximation, neglecting their finite lifetime, at a 

specific choice of impurity parameters v — 1, c — 
Solid lines are for the bands near electron-like segments of 
Fermi surface and dash-dotted lines for those near hole-like 
segments. The non-perturbed SC quasiparticle bands and 
single-impurity localized levels are shown with dashed lines. 
The narrow rectangle around the top of Sj-band delimits the 
region in Fig. 5. 



holes) of the Fermi surface, but for clarity presented in 
Fig. 3 in the same ^-reference. Besides these Sb bands, 
there appear also four (narrow) in-gap bands ±ei(±^), 
generated close to ±£o by finite concentration of impuri- 
ties, accordingly to: 

£i(0«^o + C7'|^, (24) 

As follows from Eq. 21, the Sj{^) band is located between 
its extrema Smax = + c7^£o/(A -|- Eq) at ^+ = £o + A 
and Smin = So - c7^£o/(A - eo) at = sq - A. The 
energy and momentum shifts of the extremal points by 
Eqs. 20, 21, and Fig. 3 are specific for the impurity effect 
on the multiband initial spectrum and they contrast with 
a simpler situation for an impurity level near the edge of 
a single quasiparticle band [21]. 

All these spectrum bands would contribute to the over- 
all density of states (DOS) by related quasiparticles: 
p{e) = (47rA'^)~^Im Tr Gk. The more common con- 
tributions here come from the Sb bands and they can be 
expressed through the Bardeen-Cooper-Schrieffer (BCS) 
DOS in pure crystal [26]: /9bcs(£, A) = pYsj^/e^ — A^, 



as follows: 

Pb{e) ^{\- ^f^^ P^cs (£, Ac) , (25) 

at > = A2 + 2c"f^el/{A'^ - e;^). The first factor in 
the l.h.s. of Eq. 25 describes a certain reduction of the 
BCS DOS, especially when the energy argument is close 
to the gap limits, and the shift of its gap argument is due 
to the quantum-mechanical repulsion between the band 
and impurity levels. 

More peculiar is the contribution to DOS from the Si 
bands, written as: 

Piie) « ^ , ^'-^l--^' , (26) 

at £^i„ < £2 < £^o^, and presented in Fig. 4. 

Both the effects of £& band shifts and of Si band for- 
mation can have important repercussions in the physical 
behavior of the disordered SC system and they will be 
considered below. But before this, we need to analyze 
the criteria for the considered quasiparticles to really ex- 
ist, especially in closeness to the limits of corresponding 
bands. 



V. GROUP EXPANSION AND COHERENCE 
CRITERIA 

Let us now study the crossover from band to localized 
states near the limits of Si bands, say for definiteness, its 
upper limit £max- Supposing the actual energy e < £max 
to be within the range of band states, we use the fully 
renormalized self-energy matrix, Eq. 16, up to the GE 
pair term, c^TBy, that will add a certain finite imaginary 
part Ti{$) to the dispersion law e = £i(^), Eq. 23. Then 
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FIG. 5: Parabolic approximation (dashed line) for the disper- 
sion law near the top of impurity band (solid line). 

the known loffe-Regel-Mott criterion [27, 28] for the state 
at this energy to be really band-like (also called extended) 
is written as: 



e > Vi{e). 



(27) 



To simplify calculation of the scalar function Viie)^ we 
fix the energy argument in the numerators of T-matrix 
and interaction matrices at e = EQ) obtaining their forms: 



f(e) 



■m+, i„(e) « f (e) 



£o 



(28) 



both proportional to the matrix m+ = uq® (tq + fs) 



with important multiplicative property: mA 



2rh,_|_. 



The k-summation (integration) in Eq. 28 is suitable 
done in polar coordinates over the circular segments of 
Fermi surface. Here the azimuthal integration only refers 
to the phase of numerator, resulting in a zeroth order 
Bessel function: J^" e"'=°''^de = 27rJo(a;). Since x = 
n (k-p + ^/hv-p) is typically big, x ^ 1, the asymptotical 
formula applies: Jo{x) « A/2/(7rx) cos(x — 7r/4). Then, 
for radial integration in ^ around the extremum point ^+ , 
it is convenient to decompose this fu nction in t he fast and 
slow oscillating factors: Joix) « a/2/ (-Kk^n) cos(A;+n — 
7r/4) cos[(^ — ^+) n/hv-p] with the fast wavenumber fc+ = 
kp + ^^/hvp « fcp, and to write the denominator in 
the parabolic approximation: D^{e) « (^ — ^+)^ — (5^(£), 
with 5^{e) = 4A (A + £o)^ (Emax - £) /(2c7^) (see Fig. 
5). Thus, the interaction matrix An(e) = An{e)'m+ only 
depends on the distance n between impurities, and, for e 
close to Emax, this dependence can be expressed as: 



Ar{e) 



sin kgT cos kpr, 



(29) 



where the length scales both for the monotonous decay: 

ZTT e()Pp[n.-\-e()> ^ 



27r 
kp 



c6{e) 



FIG. 6: Interaction function Ai{e) by Eq. 30 at the choice 
of parameters Smax — £ ~ 0.1 and A/ep = 5 • 10~^ dis- 
plays slow sine oscillations (solid line) and the monotonous 
envelope function (dashed line). The shadowed intervals are 
those contributing to Im B, accordingly to the condition 
{re/r)siv? keT > 1. Inset: the expansion of the rectangle 
in the main panel shows also fast oscillations by the cosine. 



and for the sine factor: k~^ = hvp/5{e), are much longer 
than /cp^ for the fast cosine. The latter fast oscilla- 
tion is specific for the interactions mediated by Fermi 
quasiparticles (like the known RKKY mechanism), unlike 
the monotonous or slowly oscillating interactions between 
impurities in semiconductors or in bosonic systems [21]. 
Now the calculation of rj(e) = c'^T{e)lmB{e) mainly 
concerns the dominant scalar part of the GE pair term: 



B{s) 



a' J a 1 



rdr 



(30) 



(since the k-dependent term in Eq. 18 turns to be negli- 
gible beside this). 

The upper integration limit in Eq. 31 corresponds to 
the condition that its integrand only has poles for r < 
Tg. In conformity with the slow and fast modes in the 
function, Eq. 30, the integration is naturally divided 
in two stages. At the first stage, integration over each 
mth period of fast cosine, around = 2iTm/kp, is done 
setting constant the slow factors, r « and sinfcer ~ 
sin fee r^, and using the explicit formula: 



Im 



dx 



1 - 4^2 ( 



Im 



(31) 



At the second stage, the summation of these results in m 
is approximated by the integration in the slow variable: 



Im 



3/2 



re sin fee 
r^/'^dr 



\/r — Te sin^ keT 



(32) 
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The numerical calculation of the latter integral results in: 



ImB = ^fik,r,), 



(33) 



where the function f{z) is zero for z < zo ~ 1.3585, and 
monotonously grows for z > Zq, rapidly approaching the 
asymptotic constant value: fas ~ 1.1478, for z ^ zg. 
Then the loffe-Regel-Mott criterion, Eq. 27, at e so close 
to Smax that k^rs '> zo, is expressed as: 



£0 a 



(34) 



and this would result in a (concentration independent) 
estimate for the range of extended states within the im- 
purity band: 



- e > Fo = 



(veo) 



3/2 



27rpF 
l + t;2 



(35) 



and its comparison with the full extension of this band, 

Smax-£min = c"/'^ {I + v'^) / {v"^ A) , would suggcst possibil- 
ity for such extended states to really exist if the impurity 
concentration surpass the characteristic (small) value: 



c » Co = 



3/2 



2v 



1 + 



(36) 



For typical values of ^ 2 eV, akp ^ 1, and A ~ 
10 meV in LaOFeAs system [5, 8, 29], and supposing a 
plausible impurity perturbation ~ 1, we estimate cq 
8 • 10~^, manifesting important impurity effects already 
at their very low content. 

However, the r. h. s. of Eq. 34 vanishes at k^re < zq, 
which occurs beyond the vicinity of the band top: 



e > F 



(37) 



Under the condition of Eq. 36, this vicinity is yet more 
narrow than Fq by Eq. 35, defining the true, even wider, 
range of extended states. 

Otherwise, for c ^ cq, the impurity band does not 
exist, then we analyze the energy range near the impu- 
rity level with the non-renormalized GE and write the 
approximate criterion for its convergence as c\B°\ <^ 1. 
This calculation is done in a similar way as before but 
replacing the interaction function, Eq. 29, by its non- 
renormalized version: 



A°{s) « ^jR^r cos kpr, 



(38) 



with fcp-Re = 27r {so/\£ — eo|) a-nd fcpro = 2£f/Co- Then 
the above GE convergence criterion is assured beyond the 
following vicinity of impurity level: 



\s - £o| » Fc = Fo exp (-Cq/^/c) , 



(39) 



defining the range of its broadening due to intcr-impurity 
interactions. The DOS function for localized states can 









max 




/ p Extended 




" states 


/Localized '•/ 





FIG. 7: Structure of the energy spectrum near the impurity 
level in function of impurity concentration. 



be only estimated by the order of magnitude within this 
range, but outside is given by: 



Pioc(e) 
Plods) 



Co^^|£-£o| 



for Fc < |£ - £o| < Fo, 



for Fo-C|£-£o|. 



(40) 



Notably, the total number of states near the impurity 
level is J pioc{s)de ^ c, alike that of extended states in 
the impurity band by Eq. 26. The summary of evolu- 
tion of this area of quasiparticle spectrum in function of 
impurity concentration is shown in Fig. 7. 



VI. IMPURITY EFFECTS ON SC 
CHARACTERISTICS 

The above results on the quasiparticle spectrum in the 
disordered SC system can be immediately used for calcu- 
lation of impurity effects on its observable characteristics. 

Thus the fundamental SC order parameter A is esti- 
mated from the modified gap equation: 



A" 



Jo 



p{e)de, 



(41) 



where A = p-pVsc is the (small) dimensionless SC pairing 
constant and the Debye energy £d restricts the energy 
range of its action. In absence of impurities, c = 0, using 
the BCS DOS in this equation leads straightforwardly to 
the known result for its non-perturbed value Ao: A~^ = 
arcsinh (£d/Ao) and thus to Aq w £De~"^/^. 

For finite c, the total DOS is combined from the contri- 
butions by the shifted main band, pb, Eq. 25, and by the 
impurity band (or level) pi (or pioc), Eqs. 26 (or 40). The 
latter contribution is ^ c, accordingly to the previous dis- 
cussion, defining a small correction beside A~^ ^ 1. But 
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a much stronger c-dependent correction comes from the 
modified main band: 



Pb{e)de « arcsinh— 



C7' 



de 



For £d Ac, the last integral is well approximated by: 
de 



C7 



/a, (£-£o)'v/A2372 A2 

with the function 

V^^ — 1 + z arccos(— 1/z) 



F{z) = z- 



3/2 



Though this F diverges at 2; 1, but actually its argu- 
ment 

Ac/^o = a/I + (1 + c/ci) with Ci = TrpF^/t;, 

is always above unity. Neglecting the small pi contribu- 
tion in Eq. 41 and taking account of the BCS relation 
A~-^ = arcsinh (ed/Ao), we express the gap equation as: 



. , Ac-Ao 
arcsmh 



cv 



Ao ci(l + w2) 



F(Ae/eo). (42) 



Its approximate solution for c ^ Ci, together with the 
relation: Ac/A = 1 + c/ [ci (l + i;^)] , lead to the desired 
expression for the perturbed SC order parameter A: 



A c l + i)2F [a/TT^ (1 + c/ci)] 

Ao ~ ci 1 + 1;2 



(43) 



that rapidly decays with impurity concentration and 
would vanish at 



c = Ci 



1 + v'^F + ^ (1 + c/ci)] 



The latter equality defines in fact a certain equation for 
c and its solution, e.g., for the above choice of u = 1, 
is c PS 0.5ci « 6 • 10~^. However, such concentrations 
would already correspond to the impurity band as wide 
as the gap itself; this goes beyond the validity of the 
above derivation and needs a special treatment (to be 
done elsewhere). 

To study another important dependence, that of the 
SC transition temperature on concentration c, one has, 
strictly speaking, to extend the above GF techniques for 
finite temperatures, but a very simple estimate can be 
done, supposing that the BCS relation A/Tc ~ 1.76 still 
holds in the presence of impurities. Then the r.h.s. of 
Eq. 43 would also describe the decay of Tc/Tcq. 

It is of interest to compare the present results with 
the known Abrikosov-Gor'kov solution for BCS SC with 
paramagnetic impurities in the Born approximation [30, 



ln(c/c - 1) 












N ^ 
N 





FIG. 8: Temperature behavior of specific heat for a SC with 
impurities presents a crossover from /3A exponent (dashed 
line) to /3eo at low enough temperature (high enough (5 = 
l/ksT). 



31]. In that approximation, the only perturbation pa- 
rameter is the (constant) quasiparticle lifetime r. In 
our framework, the can be related to lm.'E{e) at a 
proper choice of energy, e ~ [A — e| ^ A. Then, in the 
self-consistent T-matrix approximation [24] , we estimate 
~ cA/ci which leads to the relation tTc ~ ci/c, 
reaching at c > ci a qualitative agreement with the 
Abrikosov-Gor'kov universal criterion for complete SC 
suppression rTc < 0.567 (though in our case this crite- 
rion is not universal and depends yet on the perturbation 
parameter v). 

Finally, a similar analysis can be applied for the impu- 
rity effect on the electronic specific heat in the SC state, 
whose dependence on inverse temperature /3 = l/ksT is 
represented as: 



d r 



p{s)de 



(44) 



and naturally divided in two characteristic contributions, 
C = Ci + Cb, from pi and pb states: 



Ci{p) « kBC 



f3eQ 



2 cosh (/3eo/2) 



and 



CbW) « fcB(ci - c)v (/3Ac)'/' exp (-/3Ac) . 



The resulting function C(/3) deviates from the known 
low temperature behavior Co(/3) ^ exp(— /3A) for non- 
perturbed SC system at /3 > ln(ci/c — 1)/(A — eq), 
where the characteristic exponent is changed to a slower 
~ exp(— /3eo) as seen in Fig. 8. 

The same approach can be used for calculation of other 
observable characteristics for SC state under impurity 
effect, such as, e.g., differential conductivity for scanning 
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tunneling spectroscopy or absorption coefficient for far 
infrared radiation, though these issues are beyond the 
scope of this work. 

VII. CONCLUSIONS 

Resuming, the Green function analysis of quasiparticle 
spectra in an SC ferropnictide with impurities permits to 
describe formation of impurity localized levels within SC 
gap and, with growing impurity concentration, of specific 



band of extended quasiparticle states, mainly supported 
by impurity centers. Explicit dispersion laws and densi- 
ties of states are obtained for modified main bands and 
impurity bands. Further specification of the nature of all 
the states in diff'erent energy ranges within the SC gap is 
attained with analysis of group expansions for self-energy 
matrix, resulting in criteria for crossovers between local- 
ized and extended states. The developed spectral charac- 
teristics are applied for description of observable impurity 
effects. 
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